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Abstract 

In a pilot-wave theory, an individual closed system is described by a wavefunction 
ip{q) and configuration q. The evolution of the wavefunction and configuration 
are respectively determined by the Schrodinger and guidance equations. The 
guidance equation states that the velocity field for the configuration is given by 
the quantum current divided by the density | '0(g) p. We present the currents and 
associated guidance equations for any Hamiltonian given by a differential operator. 
These are derived directly from the Schrodinger equation, and also as Noether 
currents arising from a global phase symmetry associated with the wavefunction 
in configuration space. 



1 Introduction 

In the pilot-wave theory of deBroglie and Bohm [1-7], an individual closed system of 
nonrelativistic particles is described by a wavefunction ipl'Xi, . . . ,XAr,t), which satisfies 
the nonrelativistic Schrodinger equation 

( ^ \ 



^Present address. Currently Postdoctoral Fellow FWO. 
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and by particle positions Xi, . . . , xjv, for which the possible trajectories are solutions 
of the guidance equations 

where i]j = \ exp(iS'/^). Here, ■?/' is regarded as an objective physical field in configu- 
ration space, guiding the motion of an individual system. 

The 'density' |'^|^ is preserved by the particle flow, in the sense that if, over an en- 
semble of systems with the same wavefunction -0, the configurations have the 'quantum 
equilibrium' distribution . . . , Xjv, io)P time t^, then they will have the distribu- 

tion |^(xi, . . . ,Xjv,t)P at any time t. This is a consequence of the continuity equation 
for iV'p, 

WP + Ev.-(^|V.p)=o, (3) 

k=l ^ ^ 

which itself follows from the Schrodinger equation.^ 

DeBroglie arrived at the guidance equation (2), in the 1920s, by combining the 
variational principles of Maupertuis and Fermat. Bohm, on the other hand, who redis- 
covered de Broghe's theory in the early 1950s, proposed an acceleration-based equation 
of motion for the particles (Newton's equation of motion with a 'quantum potential'), 
instead of the velocity-based guidance equation. In Bohm's view, the guidance equation 
is regarded as a mere constraint on the initial momenta, a constraint that can in prin- 
ciple be relaxed [3]. Whereas, in deBroglie's view, the guidance equation is regarded as 
the fundamental law of motion. (For a full discussion, see ref. [2].) 

However, de Broghe's velocity field as given by the right-hand side of (2) cannot be 
correct in general. For while the resulting flow preserves {ipl'^ for standard Hamiltonians 
of the form appearing in (1), it will not preserve [■(/'p for Hamiltonians of an arbitrary 
form. For more general Hamiltonians, an acceptable guidance equation may be derived 
from the requirement that it preserves the density (a method that is in fact often 
used), since this property plays a key role in showing that in equilibrium the theory 
reproduces the standard quantum predictions. Given a continuity equation for a density 
It/'P, one may deflne an associated 'current'. The velocity field is then postulated to equal 
the current divided by the density. 

In this paper we derive currents, from which the guidance equations can be con- 
structed, for any Schrodinger equation such that the Hamiltonian is given by a differ- 
ential operator. More explicitly, we consider Schrodinger equations^ 

idMq,t)=H{q,V,t)ij{q,t), (4) 

where q = {qi, . . . , q^) G and V is the A^-dimensional gradient, and where H{q, V, t) 
is a differential operator, that is, a multivariate polynomial in the symbols dg-^, . . . , dqj^ 



^Note, however, that in principle the theory allows 'nonequilibrium' ensemble distributions, that is, 
distributions that differ from [8-11]. 
^Here and below we take h=l. 
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with coefficients being complex- valued functions of q and t. From (4), we derive a 
continuity equation for the density t)P with a current j'^{q,t), that is, 

dMq,t)\' + W-f{q,t)^0. (5) 

A possible guidance equation is then given by 

dq 



dt 



(6) 



Of course, one can always add a divergence-free vector field to j'^{q, t) without losing the 
validity of (5). However, it is standard to assume that the current vanish for |g| +oo. 

Note that, in the special case of one dimension {N = 1), the requirement that the 
current vanishes at infinity uniquely determines the current to be 

fM^- f dq'd^\i^{q\t)\\ (7) 

J —oo 

Other possible currents differ only by an additive g- independent term and hence do not 
vanish for |g| +oo. 

As we shall see, for a Hamiltonian that is a differential operator, we can always 
find a current j^{q) that is a local functional of the wavefunction il}{q) (that is, a j'^{q) 
that can be written as a function of q, ip{q) and ^-derivatives of ip{q) up to some finite 
order). In the uniquely-determined one-dimensional case (A^ = 1), this implies that the 
expression on the right-hand side of (7) is a local functional. On the other hand, in 
the multi-dimensional case (A^ ^ 2) one can also find currents satisfying the continuity 
equation that are not local functionals. Here is an example. Prom the Schrodinger 
equation (4) it follows that 

dtmq,t)\' + l^q,t)=0, (8) 

where 

I^{q, t) = 2Re {ir{q, t)H{q, d„ t)^{q, i)) . (9) 

We can then write 

I^{q,t) = {V-'l^{q,t)) = V • (V (V-2/^(x,t))) , (10) 

where V~^/^(g, i) = / d^ q'Gn{q—q')I'^{q' ,t) (assuming that the integral is well-defined), 
with 

{ ^\og\ql N = 2 

GM = { m/i-i) ^ > . (11) 

[ 4,r^/2|g|iV-2, iV O 

a Green's function (a fundamental solution) for Laplace's equation in A^ dimensions, 
that is, V^Giv(g) = S{q), see for example refs. [12, pp. 49-53] or [13, pp. 21-25]. The 
current 

f{q,t)^V{V-'l^{q,t)) (12) 
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then satisfies the continuity equation. This method of obtaining a current was first 
suggested by Epstein [14] (for = 3). Aithough a iocai current is arguably preferable 
when the Hamiltonian is a differential operator, nonlocal currents of the type (12) could 
be useful for other Hamiltonians, such as the positive-energy Klein-Gordon Hamiltonian 
V— + m? (which is defined by transforming to Fourier space). 

Different choices of current lead to different guidance equations and hence to different 
pilot-wave models.^ Even requiring the current to be a local functional (insofar as 
the Hamiltonian admits that) does not resolve the ambiguity, as again the current is 
determined by the continuity equation only up to a divergence- free part. Deotto and 
Ghirardi [15] considered the latter ambiguity in detail, and showed that even imposing 
standard space-time symmetries still does not uniquely fix the current. 

Before considering the general case of a multi-dimensional configuration space, we 
shall first consider the one-dimensional case. This allows us to illustrate the techniques 
involved. 



2 One-dimensional configuration space 

Consider a Schrodinger equation 

idti/j{q,t)=H{q,d„t)ij{q,t), (13) 
where g e M and where H{q,dq,t) is a differential operator, that is 

H{q, d„ t)^{q, t) = J2 hn{q, t)dy{q, t), (14) 



n>0 



where only a finite number of the complex- valued functions /i„ are non-zero. We assume 
that H is Hermitian, that is, for wavefunctions ipi and ■02, 

J dqr^H^P^ = j dq{Hij,yij2. (15) 

From (14), Hermiticity of H implies 

H{q, d„ t)ij{q, t) = 5](-9,)" {K{q, t)i;{q, t)) (16) 

(assuming that terms of the form d^{ilJihn)d^~"^~^ilJ2 vanish at infinity). As shown 
in appendix A, it then follows that the Hamiltonian is Hermitian if and only if the 
coefficients /i„ satisfy 



"^Note, however, that merely postulating a velocity field by dividing a current by the density 
might not be sufficient to obtain a satisfactory pilot- wave model (that is, a pilot- wave model that 
reproduces the standard quantum predictions). 
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From (13) it follows that 

dt\i^iq,t)\' + l^iq,t)^0, (18) 

where 

/^(g, t) = 2Re (i^*(g, t)H{q, d,, t)^l^{q, t)) . (19) 

The goal is to write I^{q,t) in the form dqj^{q,t). 

In section 2.1 we derive the current by direct calculation. In section 2.2 we derive 
it as a Noether current. We will make use of the following identity (the proof is given 
in appendix A): 

Identity 1 For any two wave functions 4>{q) and x{q)) <^'>T'd n ^ 0;^ 

cf>d^X - i-irxd^^ = d, ( I ■ (20) 

\m=0 / 

2.1 Deriving the current from the Schrodinger equation 

We have 

= i (^^*Hip - (^Hip^ 

= i^ {rhnd^ - i-iri^d^ih^r)) , (21) 

where for Hijj and (^Hijjj we have used respectively (14) and (16). Making use of 
Identity 1 (equation (20)) for = ip*hn and x = we can write I'^ — dqj'^ , with 

n-l 

/ = ij^ $^(-i)'"9-(V'*Ma,"— V- (22) 

n^l m=0 

We can also write the current in the form 

f = J2 Jnmd^^d:^r, (23) 



where the coefficients 



Jnm = i Yl (-l)^+"+'(^ ^ ^y^-^-^-'hr, n,m^O. 



(24) 



^Here and below it is understood that a sum is zero when the upper bound is smaller than the lower 
bound. 
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are completely determined by the Hamiltonian (the proof of (24) is given in appendix 
A). 

The above current is real. One way to see this is by noting that the current on 
the right-hand side of (22) (which vanishes for \q\ — > +00) equals that on the right-hand 
side of (7), which is obviously real. For the case of a multi-dimensional configuration 
space, we do not know of any simple argument that shows the reality of the obtained 
current. Instead, in appendix B we prove the reality by direct calculation, using the 
Hermiticity of the Hamiltonian. Since our discussion of the one-dimensional case serves 
as preparation for the multi-dimensional case, in appendix A we give a similar reality 
proof for the one-dimensional case. 



2.2 Deriving the current as a Noether current 

In the previous section we derived the current j'^ by direct calculation from the Schrodinger 
equation. In this section the same current is derived as a Noether current. 

As is well known, Noether's theorem states that, for theories that can be derived 
from an action principle, each continuous symmetry implies a conserved quantity. Here 
the conserved quantity is quantum probability, and the conservation law is expressed in 
terms of the continuity equation. The Schrodinger equation (13) can be derived from an 
action, and the invariance of the action under global phase transformations ijj — > e^'^xjj 
implies the conservation of quantum probability. 

Specifically, the Schrodinger equation (13) can be derived from an action principle 
with Lagrangian 

L[^, ^, r, r, t]^ Jdq (^^{r{q, t) - r{q, t)i;{q, t)) - ^iq, t)H{q, 9„ t)i;{q, t) 

(25) 

(a functional of t/j^ip^t/j* ^ip*), where ip — dt'4>- The Euler-Lagrange equations are given 

by 

5L 5L 5L 6L 

dt-j - 77 = 0, dt-j- -7-7 = 26) 

Sip Si/j S^j* Sip* 

and just yield the Schrodinger equation and its complex conjugate. 

The corresponding Lagrangian density may be taken as the function (using the Her- 
miticity of H) 

C{q, iP{q, t),r{q, t),i^{q, t), riq, t),..., t), dy*{q, t),...,t) 

= I i^riq, t)ip{q, t) - r{q, t)H{q, d^, t)iP{q, t) + c.c.) (27) 



and in terms of C the Euler-Lagrange equations (26) read 

- y (-i)na- = 0, - y{-iyd^^^- = o. (28) 
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Since the field-theoretical Noether theorem is usually discussed only for Lagrangian 
densities that depend on up to first-order derivatives of the field, we repeat the analysis in 
full. The generalization to Lagrangian densities that depend on higher-order derivatives 
is in fact straightforward. 

If L is invariant under an infinitesimal symmetry transformation ijj ^ ijj + 5ijj, 
■0* — > -0* -I- (5-0*, we have 5C — 0. We can write 5C as 

SC = -^5i> + y ——6{dy) + c.c. 



dC 
dd^ijj 



S^)+c.c., (29) 



where in the last line we have used the Euler-Lagrange equations. Using Identity 1 
(equation (20)) in the second term on the right hand side, with = dC/dd^ip and 
X = (5^, we obtain 



dC 



n-l 



and hence we obtain a conserved '2-current', with density 



in— m— 1 



dC 



= -^di^ + c.c. 



and 'spatial' current 



(5^) + c.c. (30) 



(31) 



n-l 



n^l m=0 



dC 

dd^ijj 



(32) 



(where {p^,j'^) is determined up to an additive term with vanishing 2-divergence) . 

The current associated with the global phase symmetry t/j — > t/j — iet/j, ■0* — > 
is then given by 



f = 2Re 



(33) 
(34) 



Since 



dd^i; 2^ dd^i; ' 



(35) 



7 



for n ^ 1, the current (34) can also be written as 
/ = Re 

Using d{Htjj)/dd'^'ijj = the expression in square brackets agrees with the right-hand 
side of (22), up to the constant factor e. It therefore also follows that taking the real 
part in equations (34) and (36) is redundant, as the argument of Re in those expressions 
is already real. 



n— 1 



n^l m=0 



(36) 



3 A/^-dimensional configuration space 

We now turn to the case of an A^-dimcnsional configuration space, with points q = 
{qi, ...,qN) e M^. We shall use multi-indices, which are elements of Nq , that is, we 
use an index n — {rii, . . . , tin) E N^. We define Cj G N^, for i = 1, . . . , A?", as (ej)j = 5ij 
{j — 1, . . . , N). We also define e as the element = (0, . . . , 0). We further assume 
the following standard definitions: 

n±n' — {ni± n'l, . . . , njv ± n'^^), 
n ^n' if rii ^ n'j^, i — 1, . . . , N, 
n < n' if n ^ n' and n ^ n', 





N 


\n\ - 






i=l 




N 


n\ = 


11(^^0, 




i=l 


= 






qi ' ' ' (IN 




1^ (:::)■■ 


\n'j 





where n^n' E N^. 

Consider a Schrodinger equation 



(37) 



idMq,t)^H{q,D,t)ilj{q,t), (38) 
where H{q,D,t) is a differential operator, that is 

H{q, D, t)^(g, t) = J2 hniq, t)D^'i'iq, t), (39) 



where only a finite number of the complex- valued functions /i„ are non-zero. Prom (39), 
Hermiticity of H now implies 

H{q, D, t)^{q, t) = J^i'lpD^' {K{q, t)iP{q, t)) (40) 

n>0 
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(assuming that terms of the form D^{'il)ihn)D^~^~'^^il)2-, for wavefunctions ipi and il)2, 
vanish at infinity). As shown in appendix B, it then follows that the Hamiltonian is 
Hermitian if and only if the coefficients hn satisfy 

^n-E(-i)''"'(r)^"~"^^- ^^^^ 

Prom (38) it follows that 

a*|V^(g,i)|2 + /'^(g,i)=0, (42) 

where 

7^(g, t) = 2Re (ir{q, t)H{q, D, t)^{q, t)') . (43) 

The goal is to write I^{q, t) in the form Xlili -D'^'jf (?, t). 

In section 3.1 we derive a current j^{q) by direct calculation. In section 3.2 we derive 
it as a Noether current. We will make use of the following identity (the proof is given 
in appendix B): 

Identity 2 For any two wave functions (j){q) and x(?); and n ^ 0; 
0L'"X-(-l)'"'x^> 

^-^ \ ^-^ ^ ' n l ml {■n—m—P.W /v, i v / 



1=1 \0^m^n— e. 



3.1 Deriving the current from the Schrodinger equation 

We have 



= i {rhuD^^ - i-lpi^D^iKr)) , (45) 



n>0 



where for Hip and (^Hijjj we have used respectively (39) and (40). Making use of 
Identity 2 (equation (44)) for = i>*hn and x = we can write = YliLi D^^jf, with 



E (-l)H^^ - ''% D-{rK)D----^^ip. (46) 

(up to a divergence- free term). In appendix B we prove that this current is real. 
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We can also write the current in the form 

jf = E Ji,nmD''^|^D"^r, (47) 

where the coefficients 

with n,m ^ 0, are completely determined by the Hamiltonian (the proof of (48) is given 
in appendix B). 



3.2 Deriving the current as a Noether current 

The current can also be derived as a Noether current associated with the global phase 
symmetry of the Lagrangian 



where ip = dtip. The Euler-Lagrange equations are given by 

^6L 5L ^ 6L 5L 

Sip oi/j Sip* Sip* 

which yield the Schrodinger equation (38) and its complex conjugate. 
The corresponding Lagrangian density may be taken as the function 

£ = i {i{ip*i) - ip*ip) - iP*HiP - (HiP) * ip) . (51) 



2 

In terms of £, the Euler-Lagrange equations (50) read 

d^tP ' dD^iP ' 'dip* ^ dD^iP* ^ ' 

If C is invariant under an infinitesimal symmetry transformation ip ^ ip + Sip, 
ip* ^ ip* + Sip*, we have S£ = 0. We can write S£ as 

SC = ^Si, + y ^^S(D"i-) + CO. 

- (I**) + E {-^^'"m - (-1)""^" {-^) + (63) 
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where in the last hne we have used the Euler-Lagrange equations. Using Identity 2 
(equation (44)) in the second term on the right hand side, with = d£/dD"-ijj and 
X = Sip , we obtain 



= 6£ 



X 



n\ 


\m\ 


\n 


— m 


-e,|! 


\n\\ 


m\ 


in 


— m 


-e,)! 



\m\ 



dC 



D"-^'5iP + c.c. (54) 



and hence we obtain a conserved '(A^ + l)-current' (up to an arbitrary additive term 
with vanishing (A^ + l)-divergence), with density 



p'f' = -jSi/j + c.c. 
dip 



(55) 



and 'spatial' currents 

A* = E E (-1) 



n! 


m ! 


\n 


— m — 


e.|! 


n ! 


m\ 


[n 


— m — 


ei)\ 



D"-^'8^^ c.c. (56) 



The current associated with the global phase symmetry ip ip — ieip, ip* — > il)*+ieip* 
is then given by 



p^^e\iP\', 



jf = 2Re 



Since 



(57) 



--EE (-i: 



|m| 





m 


I 


n — m — Ci 


I 


\n 


I 


m 


(n — m — Cj)! 



^ip 



dC _ I djHiP) 
dD^ip ~ 2^ dD'^ip ' 

for n > 0, the current (58) can also be written as 



(58) 
(59) 



Re 



-EE (-1)""' 



X 



n! 


m 


I 


n — /;/ — Cj 


I 




I 


m 


(n — m — Cj)! 



D*" 1p 



dD 



'"^il) j 



D 



n—m—ei 



iP 



(60) 



Using d{Hip)/dD^ip = hn, the expression in square brackets agrees with the right-hand 
side of (46), up to the constant factor e. It therefore also follows that taking the real part 
in equations (58) and (60) is again redundant as the argument of Re in those expressions 
is already real. 
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4 Comparison with related work 



The construction of quantum probability currents and associated guidance equations, 
for certain classes of operators, has been considered before by a number of authors, in 
particular Brown and Hiley [16,17], Diirr et al. [18,19], and Gambetta and Wiseman 
[20] .6 

4.1 Brown and Hiley 

Brown and Hiley consider one-dimensional Hamiltonians of the form 

H^J2^-^nq"T, (61) 

m,n 

where the h^n are constant and p is the momentum operator conjugate to q. In the 
configuration representation, where {q\q\q') = qS{q — q'), {q\p\q') = —idgS{q — q'), the 
Hamiltonian operator is given by 

{q\H\q') = H{q, d,)6{q - g') = E hmnq"'{-id,rS{q - q'), (62) 

m,n 

which corresponds to a Hamiltonian of the form (14) with /i„ = ^mn?"^(~i)"- 
They find a current that can be written in the compact form 

f{q,t)^{q\dp{pm)\q), (63) 

where p{t) = \ip{t)){ip{t)\ is the density operator, and dp is the symbolic differential 
operator of Born and Jordan [16, 17, 22], which acts on pq^'p^ as 

n 

dpipcTP^) = (64) 

and whose action is extended to pH by linearity. It can easily be shown that, in its 
domain of validity, this current coincides with that given in (22). 

It is straightforward to verify that Brown and Hiley's expression for the current also 
applies to Hamiltonians of the form 

H = J29n{q,t)r, (65) 

n 

^Stone [21] considers Hamiltonians of the form (14) where the /i„ are constant. Instead of writing 
the guidance equation in terms of derivatives of ijj and ip* , Stone attempts to write it in terms of the 
real and imaginary parts of {idq)'""ilj / ip . However, Stone's expression for the guidance equation appears 
to be incorrect. While it is rather straightforward to rewrite the guidance equation in terms of Stone's 
variables, this does not seem particularly useful. Therefore we do not present it here. 
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where (g|^„(g, t)\q') = Qniq, t)6{q — q'), and which correspond to Hamiltonians H{q, dq, t) 
of the form given in (14) with /i„ = i)". (The Born- Jordan derivative acts on 
P9n{Q, t)p" in the same way as on pq^p^, replacing by 'gn{q, t) in (64).) 

While we have not investigated this further, a similar form can probably be obtained 
in the case of a multi-dimensional configuration space. 

4.2 Diirr et al., Gambetta and Wiseman 

Diirr et al. [18, 19] and, independently, Gambetta and Wiseman [20] consider multi- 
dimensional Hamiltonians of up to second order in the momenta, that is Hamiltonians 
of the form (we follow a notation similar to that of Gambetta and Wiseman) 

H aijPiPj + ^ biPi + c, (66) 

where aij,bi,c arc arbitrary time-dependent functions of the g^. In the configuration 
representation, these correspond to Hamiltonians of the form (39) with the /i„ zero for 
|n| > 2. 

They find a current that can be written as 

jt{q,t) = Rei{m\l){Qm)\m)) = Re(g|t;.(t)p(t)|g), (67) 

where % — i[H, qk] is the velocity operator. This expression for the current is not 
valid for Hamiltonians containing higher-order terms in the momentum operators, cor- 
responding to some non-zero /i„ for \n\ > 2. However, in its domain of validity, it can 
be easily be shown to coincide with the current given in (46). 

5 Conclusion 

We have presented de Broglie-Bohm guidance equations for any Schrodinger equation 
such that the Hamiltonian is given by a differential operator. These results may be ap- 
plied to the development of possible pilot-wave interpretations for quantum theories with 
a Schrodinger time evolution. While we have considered only finite-dimensional config- 
uration spaces, the equations carry over — at least formally — to field-configuration 
spaces. Thus the work can also be applied to the development of pilot-wave interpreta- 
tions for quantum field theories, in terms of field ontologies. 
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A Proofs for a one-dimensional configuration space 

Proof of the Hermiticity condition (17). Prom (16), and using Leibniz' rule, we have 



H{q,d„t)i^iq,t) = ^f;(-i)"('^')a;^-/i:(g,t)a-v(g,^)- 

n^O m=0 ^ ^ 

Using the identity ^^^^ Y.l=o fnm = En>o Em>n fmn we have 



(68) 



H{q,d„t)i^{q,t) = J2 E(-ir(^)^r"^m(^'^)^.X^'^)- (69) 



Comparing this expression with (14) gives 



(70) 



Conversely, the condition (70) imphes that 

J2 hn{q, t)d:^^{q, t) = 5^(-a,)" {K{q, t)^{q, t)) (71) 

SO that the Hamiltonian is Hermitian. □ 

Proof of Identity 1. The identity is proved by working out the right-hand side of equation 
(20): 

(n— 1 \ n—1 

Y,{-ird^(t>d-^-^-\ = Y,{-ir (srv^r""'^ + ^rH"""^) • (72) 
m=0 / m=0 

Using the identities 

n-1 n-2 

^(-1)™5™+V9r"^-'X = E(-l)™5;"+V5r'""'x - {-lTxd;ct> (73) 

m=0 m=0 

and 

n—1 n—2 

^(-i)-9-09^"— x= E (-i)"^^'9rv9r""'x 

m=0 m=— 1 

n-2 

= - J](-i)"^a™+V5r™"'x + 05;%, (74) 

m=0 
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we find that (72) becomes (pd^x ~ {~^)"'x9q(p, which is the left-hand side of equation 
(20). □ 

Proof of equation (24). Regarding as a function of the variables d^ipd^ip*, we have 
that Jnm = dj'l' /dldgipd^ip*). To calculate dj"^ /d{dgiljd^ilj*), we first use Leibniz' rule 
to write 

r-l 



= iEEE(-i)^f !)^5^*^rMr"V- (75) 

r>l s=0 t=0 ^ ^ 



We then have 



df 



did:^^d-r) = ^ E E t(-^r (j) ^^%s.,r-s-^Sr.,t 

\ q r J \ / 

= i E E(-l)^(^)C"^Mn,.-.-l 



=i E {-ir''^'("~ll~^)dr-"'-'hr. (76) 



□ 



Proof that the current ^^zven 6y (22) is rea/. The proof is best given by starting 
from the form (23) of the current. The current is real if and only if = J^^ (for all 
n,m ^ 0). 

Before proving the latter, we first show that: 



s=n+m+l 



for r ^ n + m-t- 1, with r,n,m G N. We prove this by induction on I = r — n — m — 1 ^0. 
For Z = 0, we have that both sides of equation (77) are equal to (— l)''m!/r!. Suppose 
now I > and suppose the identity holds for all smaller values of I. Using the identity 

r — n — m — 1\ fr — n — m — 2\ fr — n — m — 2\ 

= + 1 K 78 

r — s J \ r — s J \ r — s — 1 J 

which is vahd for any r,n,m and s [23, p. 174],^ the sum on the left-hand side of (77) 

'^We use the definition of binomial coefiicients in ref. [23, p. 154], which reads 



m>0 

rn(m— ' 

1, m = (79) 

0, m < 



for integer n, m. 
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decomposes into two sums, which we call Si and 82- Let us first consider 



c — -n -l-'^n 4- 1 ^ ' 



s=n+m+l 

Since ('^"'^l^"^) — ^iors — n + m+1 and r — n — m — 1>0, we have 



s\ \ r — s 

s=n+m+2 



Applying the induction hypothesis we find 



Similarly we have 



s — n — l)\fr — n — m — 2 
r — s — 1 



s=n+m+l ^ 

g ^_^^^{s-n-iy. fr-n-m-2 



s\ \ r — s — 1 

s=n+m+l 

.„+„+im!(r-m-2)! 



r!n! \ r — m — 1 r — m — 1 



rW. 



(80) 



5.= E (-ir^^^r"'::r'y 



S, = (-1^""" + ^'"';-'"-^)'. (82) 



Combining these results we find that the left-hand side of equation (77) equals 
S, + S, = ^-un^rn^^ rnlir - m - ly. ( m + 1 ^ 



(84) 



which is equal to the right-hand side of (77). This completes the proof of the identity 
(77). 

Consider now 



r'^n+m+1 ^ ' 



(85) 



Using respectively the Hermiticity condition (17), the identity '^r>n+m+i'^s>r frs — 



Z)r^n+m+iX)s=n+m+i/sr-' ^he factorial form of the binomials, and the identity (77), we 
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find 



= i E E (-1)— ^ 

r^n+m+l s=n+m+l 



\^ ^ ^ m!(r-n-m-l)! « 

X 



s=n+m+l ^ ' 



= ^^n- (86) 

□ 



B Proofs for an A^- dimensional configuration space 

Proof of the Hermiticity condition (41). From (40), and using Leibniz' rule, we have 



H{q,D,t)i;{q,t) = J2 E i'^P (l^) D''-'^h:{q,t)D^i;{q,t). 
Using the identity ^^^^ Eo^m^n /nm = En^o Em^n /mn we have 
H{q, D, t)i;{q, i) = E E ("l)'™' Q 



(87) 



(88) 



Comparing this expression with (39) gives 



J](-l)HMi^— /i* 



(89) 



Conversely, the condition (89) imphes that 

J2 hn{Q, ^)^X9, t) = Yi-lpD'^ (Kiq, t)ij{q, t)) , (90) 

so that the Hamiltonian is Hermitian. □ 
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Proof of Identity 2. The identity is proved by working out the right-hand side of (44) : 

^-^ \ ^-^ \n\\ ml (n — m — ej! / 

^ I I I! I __ 11 

E E (-l)'"!^^ y (D"-+-'0D"-"--^x + D^<j>D"x) . (91) 



i=l O^m^n— e, i i v / 

We can write 

m! (n — m — Cj)! 

^-^ (m — ej! In — mj! 

Iml! In — ml! m. 



y (_l)|m|+i K^^L - Ll^D^ct>D"x 
^-^ ml {n — my. \m' 

m\\ \n — m\\ m^ 
m\ {n — m)\ \m\ 



^ ^ m\ in — m)\\m 

' ^ Tn\ (n — mil Iml 

J2 (-l)H+i^ _ {-lp\^^D-<j>x. (92) 

ml (n — m)l \m\ nl \n\ 



0<m<n 



(We can extend the range of the sum in the second-last hne, because the terms with 
mi — are zero.) Similarly, we can write 



^-^ ml [n — m — Cj)! 



m ! \n — m — e,- 
^-^ ml (n — m)l \n — m\ 

, , Im.ll In, — mil (n.^ — mA 



0<m<n 



m ! 


n 


— m ! 




- mi) 


m\ 


(n 


— m)\ 


n 


— m\ 



m ! 


n 


— m\\ 


(Ui 


- rrii) 


ml 


(n 


— my. 


\n 


— m\ 



D'^(t)D"x 

D'^(t)D"x + ^ f^'/'^^X- (93) 



\n\ 


! Hi 


n\ 


\n\ 



0<m<n 

Using (91), (92) and (93), the right-hand side of (44) becomes 

<AD»x - (-l)l"ix£>> 

+ E E (-l)'""o^j^(^;^-f^)i''»*^"-X. (94) 

* ^ i ^ In 1 1 nm \ in — m 11 \ In — m I \Tfl\ ' 



i=l 0<m<n 
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Using 



y(^-^^^^-^)=0, (95) 
\ \n — ml \m\ ' 



1=1 

(94) is equal to (j)D''x - (-l)'"'x-D"0, which is just the left-hand side of (44). □ 

Proof of equation (48). Regarding the jf as functions of the variables D'^ipD'^ip* , we 
have that Ji^nm = djf /d{D''ijjD'^i;*). To calculate 9jf /^(D^D"*^*), we first use 
Leibniz' rule to write 



= 'E E Ef-D'-'^^j^^^^HrO^v^'-''-^'-'-'^- P") 



We then have 

'E E E' 



diD^i^Dmr) h^jy_ £f2 ^ \r\\ s\ {r - s - ei)\tr " 



r—s—ei'Jm.t 



r !(r — n — ej! n! V / 

(97) 
□ 

Proof that the currents jf given by (46) are real. The proof is best given by starting 
from the form (47) of the currents. The currents are real if and only if Ji^nm — Ji,mn i^^^ 
all n, m ^ 0). 

Before proving the latter, we first show that: 

E (-^)'"M^C"7-r 1 - '-^)""""" ''""'h.n!""" - (^«' 

n+m+ei-^s-^r I I \ / I I I I 

for r ^ n+m+Ci, with r,n,m G N^. We prove this by induction on / = r—n—m—ei ^ 0. 
For Z = 0, we have that both sides of equation (98) are equal to (— l)l'"l|m|!/|r|!. Suppose 
now I > and suppose the identity holds for all smaller values of I. Since I > there 
exists a. j e {1, . . . , N} such that I — Sj— r — n — m — Si — Sj'^O. Using the identity 

fr — n — m — eA /r — n — m — ei — eA /r — n — m — ei — e A , , 

H + ' 99 

\ r — s J \ r — s J \ r — s — ej ) 
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which is vahd for any r,n,m and s [23, p. 174], the sum on the left-hand side of (98) 
decomposes into two sums, which we call 5*1 and 5*2. Let us first consider 



E (-i)'-' ""!;"'^' r (100) 

Since = for Sj = rij + rrij + Sij and r — n — m — Ci — Cj ^ 0, we have 



n+m+Si+ej ^s^r 

Applying the induction hypothesis we find 



r ! n ! 



Similarly we have 



s — n — ei\ \ / r — n ~ m — Ci — Cj 



r — s — Cj 

' s — n — ei\ \ I r — n — m — e-i — 
r — s — Cj 



= (_l)|n+H+il^n;-^-^--^^-|'. (103) 
|r — ej|!|n|! 

Combining these results we find that the left-hand side of equation (98) equals 

|r|!|n|! \ |r — m — ej| |r — m — e 

^ (_i)in+H+iM^_:_^_zM, (104) 

which is equal to the right-hand side of (98). This completes the proof of the identity 
(98). 

Consider now 

J..n^ = i E r " " " '1' ^ - " - D'-'—"K (105) 



Using respectively the Hermiticity condition (41), the identity '^r^ri+m+ei'^s>r f^s = 
Hr^n+m+ei'^n+m+ei^s^r fsr^ the factorial form of the binomials, and the identity (98), 
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we find 

r! \r — n — ej\\ \n\\ I r — n — ei \ ( s 

I "I \ I ' T^'t-r-^ I -r-L /.111! /, 1 I 

r^n+m+Ci s^r 



^ s\ \s-n-e,\\ \n\\ - n - eA /r\ 

\s\\ {s — n — ei)l nl \ m J \s J ^ 

= i y (-l)l"+"l+'-r-r7 — -.D''-"-^^ 

^^-^ n\m\{r — n — m — ei)\ 

X ^ ( ^y,| |g-n-e^|! /r-n-m-eA 

_• ( -^\r^m\ ^- \'c — — tj^. |r?7-|! (r — m — eA 

" .^nT;;:+e. |r|!(r-m-eO! m! V ^ 7 

= (106) 

□ 
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